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Abstract 

In this paper we investigate the quantum Zeno and anti-Zeno effects without using any particular 
model of the measurement. Making a few assumptions about the measurement process we derive 
an expression for the jump probability during the measurement. From this expression the equation, 
obtained by Kofman and Kurizki [Nature (London) 405, 546 (2000)] can be derived as a special 
case. 
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I. INTRODUCTION 



The description of the measurement process has been a problem since early development 
of quantum mechanics jl|. During recent years the measurement problem attracted much 
attention due to the advancement in experimental techniques. Nevertheless, the full un- 
derstanding of quantum-mechanical measurements has not been achieved as yet. Typically, 
the measurement in quantum mechanics is described by von Neumann's state reduction 
(or projection) postulate However, this postulate refers only to an ideal measurement, 
which is instantaneous and arbitrarily accurate. Real measurements are represented by the 
projection postulate only roughly. 

The so-called "quantum Zeno effect" is directly related to the measurement problem. In 
quantum mechanics the short-time behavior of nondecay probability of an unstable particle 



is not exponential but quadratic |2|. The deviation from the exponential decay has been 
observed by Wilkinson et al. jjj. Using the behavior of nondecay probability Misra and 
Sudarshan in 1977 showed, that the frequent observations can slow down the decay. 
An unstable particle would never decay when continuously observed. Misra and Sudarshan 
have called this effect the quantum Zeno paradox or effect. The very first analysis does 
not take into account the actual mechanism of the measurement process involved, but it is 
based on an alternating sequence of unitary evolution and a collapse of the wave function. 
The quantum Zeno effect has been experimentally proved [5| in a repeatedly measured two- 
level system undergoing Rabi oscillations. The outcome of this experiment has also been 



explained without the collapse hypothesis 



Later it was realized that the repeated measurements could not only slow down the 
quantum dynamics but the quantum process may be accelerated by frequent measurements, 
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15| . This effect was called a quantum anti-Zeno effect. Quantum 
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Zeno and anti-Zeno effect were experimentally observed in an atomic tunneling process 

Simple interpretation of quantum Zeno and anti-Zeno effects was given in Ref. [ljj. Using 
projection postulate the universal formula describing both quantum Zeno and anti-Zeno 
effects was obtained. According to Ref. [ill , the decay rate is determined by the convolution 
of two functions: the measurement-induced spectral broadening and the spectrum of the 
reservoir to which the decaying state is coupled. 

In this paper we analyze the quantum Zeno and anti-Zeno effects without using any 
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particular measurement model and making only few assumptions. We obtain a more general 



expression for the jump probability during the measurement. Expression, derived in Ref. [11 1 
is a special case of our formula. 

The work is organized as follows. In Sec. [HI we present the description of the measure- 
ment. A simple case is considered in Sec. 11111 In Sec. II VI we derived a general formula for 
the probability of the jump into another level during the measurement. The pulsed measure- 
ments when there is a period of the measurement-free evolution between the measurements 
is analyzed in Sec. [V[ Particular case of the expression, obtained in Sec. IIVI is investigated 
in Sec. IVII Section fVIII summarizes our findings. 



II. DESCRIPTION OF THE MEASUREMENT 



We consider a system that consists of two parts. The first part of the system has the 
discrete energy spectrum. The Hamiltonian of this part is Hq. The other part of the system 
is represented by Hamiltonian H 1 . Hamiltonian Hi commutes with H . In a particular 
case the second part can be absent and Hi can be zero. The operator V(t) causes the 
jumps between different energy levels of H . Therefore, the full Hamiltonian of the system 
is of the form H$ = H + Hi + V(t). The example of such a system is an atom with the 
Hamiltonian Hq interacting with the electromagnetic field, represented by Hi, while the 
interaction between the atom and the field is V(i). 

We will measure in which eigenstate of the Hamiltonian H the system is. The measure- 
ment is performed by coupling the system with the detector. The full Hamiltonian of the 
system and the detector equals to 

H = H s + H D + Hi, (1) 

where Hjj is the Hamiltonian of the detector and Hi represents the interaction between the 
detector and the measured system, described by the Hamiltonian Hq. We can choose the 
basis \na) = \n) ® \a) common for the operators Hq an d Hi, 

H Q \n) = E n \n), (2) 
Hi\a) = E a \a), (3) 

where n numbers the eigenvalues of the Hamiltonian H and a represents the remaining 
quantum numbers. 
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The initial density matrix of the system is ps(0). The initial density matrix of the detector 
is Pd(0). Before the measurement the measured system and the detector are uncorrelated, 
therefore, the full density matrix of the measured system and the detector is p(0) = ps(0) <E> 
Pd(0). The duration of the measurement is r. 

When the interaction of the detector with the environment is taken into account, the 
evolution of the measured system and the detector cannot be described by a unitary op- 
erator. More general description of the evolution, allowing to include the interaction with 
the environment, can be given using the superoperators. Therefore, we will assume that 
the evolution of the measured system and the detector is given by the superoperator S(t). 
The explicit form of the superoperator S(t) can be obtained from a concrete model of the 
measurement. 

Due to the finite duration of the measurement it is impossible to realize the infinitely 
frequent measurements. The highest frequency of the measurements is achieved when the 
measurements are performed one after another without the period of the measurement- 
free evolution between two successive measurements. Therefore, we model a continuous 
measurement by the subsequent measurements of the finite duration and finite accuracy. 
After N measurements the full density matrix of the measured system and the detector is 

p(Nr)=S(rfp(0)- (4) 

We assume that the density matrix of the detector, pu(0), is the same before each mea- 
surement. Such an assumption is valid when the initial condition for the detector, modified 
by the measurement, is restored at the beginning of each measurement or each measurement 
is performed with a new detector. For example, the detector can be an atom which is ex- 
cited during the measurement. After the interaction of the atom with the measured system 
is interrupted, the atom returns to the ground state due to spontaneous emission, and the 
result of the measurement is encoded in the emitted photon. Thus the initial state of the 
detector is restored. 

III. MEASUREMENT OF THE UNPERTURBED SYSTEM 

In this section we investigate the measurement of the unperturbed system, i.e., the case 
when V(t) = 0. 
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We assume that the measurement of the unperturbed system is a quantum non-demolition 



measurement 
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The measurement of the unperturbed system does not change 
the state of the measured system when initially the system is in an eigenstate of the Hamilto- 
nian H Q . After such an assumption, the most general form of the action of the superoperator 
S(t) can be written as 

S(r)[\na)(ma'\ ® p D (0)} = \na)(ma , \e^^ T <g> S na>maf (r)p D (0), (5) 



where 



Uma'.na — ^{E m + E a i — E n — E a ) (6) 

n 



and the superoperator S na>ma i (r) acts only on the density matrix of the detector. The full 
density matrix of the detector and the measured system after the measurement is 

p{r) = <S(r)p(0) = \na)(p s )n a , ma 'e %Uma '> naT (ma'\® S na)ma i(T)p D {Q). (7) 

nct,ma' 

From Eq. (J7J) it follows that the non-diagonal matrix elements of the density matrix of the 
system after the measurement (ps)na,ma'( r ) are multiplied by the quantity 

Fn a ,rna'(r) = Tl {S na , ma , {f)p D (0) } . (8) 

Since after the measurement the non-diagonal matrix elements of the density matrix of the 
measured system should become small (they must vanish in the case of an ideal measure- 
ment), F na)ma t{r) must be also small when n ^ m. 



IV. MEASUREMENT OF THE PERTURBED SYSTEM 

The operator V(t) represents the perturbation of the unperturbed Hamiltonian Hq + H\. 
We will take into account the influence of the operator V(t) by the perturbation method, 
assuming that the strength of the interaction between the system and detector is large and 
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the duration of the measurement r is short. Similar method was used in Ref. 

We assume that the Markovian approximation is valid, i.e., the evolution of the mea- 
sured system and the detector depends only on their state at the present time. Then the 
superoperator S, describing the evolution of the measured system and the detector, obeys 
the equation 

= £{t)S, (9) 



S^(t,t )=Tex P ^C (t')dt^ (13) 



where £ is the Liouvilian. There is a small perturbation of the measured system, given by 
the operator V. We can write £ = Co + £y, where £y is a small perturbation. We expand 
the superoperator S into powers of V 

S = S<®+SW+ 5^ + -.. (10) 

Then from Eq. JOJ) it follows 

J^(°> = Co{t)S®, (11) 

J^<S« = £ (t)5» + £v(t)5^. (12) 

We will denote as S^°'(t,t ) the solution of Eq. with the initial condition S^°'(t = 
to, to) = 1. The formal solutions of Eqs. (fTTj) and (|T2"j) are 

1 1 
and 

S«(f,0) = PdtxS^tJCvitJSWfaO). (14) 

JO 

Here T represents the time-ordering. In the second-order approximation we have 

S(t,0) = 5 (0) (t,0)+ /*d*i5^(t,ti)£ v (ti)5 (0) (ti,0) 

jo 

+ /V / -tl ^ 2< S( )(t,t 1 )£ v (t 1 ) < S(°)(t 1 ,t 2 )£ y (i 2 )«S( )(t 2 ,0). (15) 
jo jo 

Using Eq. (fTUj). the full density matrix of the measured system and the detector can be 
represented as 

p(t) = pW(t)+pU(t)+pW(t) + ... , (16) 

where 

pW(t) = 5»(t,0)p(0). (17) 
Let the initial density matrix of the system and detector is 

p(0) = \ia)(ia\®p D (0). (18) 

The probability of the jump from the level \ia) into the level \fat') during the measurement 
is 

W(ia -> fa') = Ti{\fa')(fa'\p(r)}. (19) 



Using the equation (J3J) we can write 

S (0) (Mo) [\na)(ma'\®p D (0)] = \na){ma'\e^,nJ «Si°J )W (t,to)^(0). (20) 

From Eq. (|2T)j) it follows that the superoperator S^' ama with the equal indices does not 
change the trace of the density matrix /3d, since the trace of the full density matrix of the 
measured system and the detector must remain unchanged during the evolution. 

When the system is perturbed by the operator V{t) then the superoperator C v is defined 
by the equation 

C v (t)p=±r[V(t),p). (21) 

The first-order term is pW(t) = S^(t, 0)p(0). Using Eqs. (JH), (HHJ), flUJ), and (J2U), this 
term can be written as 

= Li r^fba!)^^^)^^--!^ 2 )^!®^^^^) 

- [i«)^ pcei (t 2 ) e i -— ^^(pce!! ® <Sg pcei (t,t2)) <5^(t 2 , 0)p^(0). (22) 

When i ^ f then the first-order term does not contribute to the jump probability, since from 
Eqs. (I19|) and (J22j) it follows that the expression for this contribution contains the scalar 
product (fa'\ia) = 0. 

For the second-order term p^ 2 '(t) = S^ 2 >(t, 0)p(0), using Eqs. (fTlj) and (j2~Uj) . we obtain the 
equality 

Tr{\fa'){fa'\pW(t)} = i jf* c^Tr {(/^(^P^OIM " (/a'|p (1) (ti)V^(ti)|/a'>} . 

(23) 

In Eq. f)23p the superoperator Sf2>,fa' * s om itted, since it does not change the trace. Then 
from Eqs. (|2^|) and (|2~Hj) we obtain the jump probability 

W(ia^fa') = ^ f Q dt x J^dhTr^Vfa^MViajMS^ 
+ V7 Q ^ iQ (t 2 )U iaJa K^)4°',-( t i' t 2)e^-^' (tl -* 2) ) 

xC(«».°)M0)}- ( 24 ) 

Equation (}2"^|) allows us to calculate the jump probability during the measurement when 
the evolution of the measured unperturbed system is known. The explicit form of the 
superoperator S^ ma , can be obtained from a concrete model of the measurement. The main 
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assumptions, used in the derivation of Eq. (JHJ), are Eqs. © and ©, i.e., the assumptions 
that the quantum measurement of the unperturbed system is non-demolition measurement 
and that the Markovian approximation is valid. Thus, Eq. (|24|) is quite general. 
The probability that the measured system remains in the initial state \ia) is 

W{ia) = 1 - J2 W ( ia -> /«')■ ( 25 ) 
/,"' 

After N measurements the probability that the measured system remains in the initial state 
equals to 

W{ia) N « exp(-RNr), (26) 

where R is the jump rate 

R = ^2-W{ia -> fa'). (27) 

f,a' T 

V. FREE EVOLUTION AND MEASUREMENTS 

In practice, it is impossible to perform the measurements one after another without the 
period of the measurement-free evolution between two successive measurements. Such inter- 
vals of the measurement-free evolution were also present in the experiments demonstrating 



the quantum Zeno effect j^, I16J, |22j. Therefore, it is important to consider such measure- 
ments. This problem for the definite model was investigated in Ref. j^j]. 

We have the repeated measurements separated by the free evolution of the measured 



system. For the purpose of the description of such measurements we can use Eq. 
obtained in Sec. IIVI The duration of the free evolution is t f , the duration of the free 
evolution and the measurement together is r. The superoperator of the free evolution 
without the perturbation V is S F \t), the superoperator of the measurement is S^(t,to). 
We will assume that during the measurement the superoperator Co does not depend on time 
t. Then the superoperator Sj$(t,to) depends only on the time difference t — t . Therefore, 
we will write S^\t — t ) instead of S^(t, to). When the free evolution comes first and then 
the measurement is performed, the full superoperator equals to 



c(o) (f f 

,J na,ma'\ L i L ^ 



^Mn a , m a'(t ~ h), t > h > t f and T > t > ti, 

S ( p\t - t{), t f > > and r F > t > t u (28) 

S Mna ma'it ~ t f)<SP ( T F ~h), T F > t X > and T>t>T F . 
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Equation (|28|) can be written as 

s£w(Mi) = 5£ ) nQima ,(i-t 1 )e(t 1 -r F )+4 0) (*-*i)©(rF-t) 

+<s£L,w(* - r F )SP(r F - t ± )e(t - t f )Q(t f - t,), (29) 

where B is Heaviside unit step function. From Eqs. (J24j) and (J29|) it follows that the jump 
probability consists of three terms 

W(ia fa') = W M (ia fa') + W F (ia fa') + W^ia -> fa'), (30) 

where the jump probability during the free evolution is 

W F (ia -> fa') = ^ r dt x r dt 2 V fa/> Uh)V iaJa/ (t 2 )e l ^'^ h - t2 \ (31) 
n Jo Jo 

the jump probability during the measurement 

W M {ioc^fa') = ^ [ T dtx f 1 dt 2 Ti{(v }al ^ti)V ia ja\t 2 )S { S iaJa ^ 

ft Jrp Jrp *■ ^ 

+ v fa ,,ut2)v iaJal {ti)s { S fal , a {ti -t 2 y^>^-^) 

S { M\ a , ia {t 2 -r F )Sf{r F )p D {Q)}, (32) 



x 



and the interference term is 



W T (ia -> fa') = ^ £ dh fj dhTr^Vfa^MV^i^S^ 

+ V^^V^jAh^^^h - TF )e^f^-^) sP(r F )p D (0)} .(33) 

If we assume that the free evolution does not change the density matrix of the detector 
and the perturbation V does not depend on time, we have the jump probability during the 
measurement-free evolution 

4 sin 2 ( huJf a i i a T F ) 

W F (ia - fa') = \V iaJa ,\ 2 -AiJ-J J -, (34) 

the jump probability during the measurement 

W M (ia^fa') = hv iaJa ,\ 2 T dh f 1 dhTiUsf^^h-t^^^-^ 
Ti Jtf Jtw 



xS ( 2 ia>ia (t 2 -r F )p D (0)}, (35) 



and the interference term 



W z (ia -> fa') = \Vi 12 



2 sin (hujfa' 



ia, fa' | -2 

ft UJf a ',ia 

x £ dt! Tr { (s^^fa - r F )e^'-(^-M 

+ SfiW** " ^e^'^-M) p D (0)} . (36) 

VI. SIMPLIFICATION OF THE EXPRESSION FOR THE JUMP PROBABILITY 

The expression for the jump probability during the measurement can be simplified if the 
operator V does not depend on time t. Then Eq. f!24j) can be written as 

W ( m ^ fa ') = ^\V iaJa ,\ 2 Re jT dt x J* dt^f*'^'^ Tr{s£] fa ,(ti, t 2 )S[% a {t 2 , 0)p D (0)}. 

(37) 



Introducing the function 



G(uj)f a ' >ia = \V ia j Q '\' 2 5 ( j\E a ' — E a ) — 10 ) (38) 



we can rewrite Eq. (}3~Tj) in the form 



2ttt f°° 

W(ia -> fa') = -tt G(uj) fa/i iaP{uj)i a ja'duj, (39) 

lL J — oo 

where 

P(uj) ia , fa , = — Re r* 2 e^^^-«^{5^(ti ) * 3 )5g fa (« a> 0)p D (0)}. (40) 

7TT JO JO 

Equation ()39|) is similar to that obtained by Kofman and Kurizki in Ref. 

Further simplification can be achieved when the superoperator £ does not depend on 
time t and the order of the superoperators in the expression Tr ^S^j a ,(ti, t 2 )S^ ia (t 2 )p D (0)^ 
can be changed. Under such assumptions we have 

^{S^ a/ (t u t 2 )sS M (t2)p D (0)} = Ti{slZ a (t 2 )S^ f Ati,t 2 )p D (0)} = F ia j a i{ti — t 2 ), 

(41) 

where F ia j a >(t) is defined by Eq. (jHJ). After changing the variables into u = t\ — t 2 and 
v = ti + t 2 from Eq. (f4*Uj) we obtain 

P{u)ioL,fct> = - Re / (l - -) F ia f a '(u) exp (i(a> - c^/)m) du. (42) 
7r jo V r/ 
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A. Decaying system 



We consider a decaying system with the Hamiltonian Hq that due to the interaction with 
the field decays from the level |i) into the level |/). The field initially is in the vacuum state 
\a = 0). Only the energy levels of the decaying system are measured and the detector does 
not interact with the field. Then Sf^j a , and P(uj)i a j a > do not depend on a and a'. Using 
Eqs. fl2Tj) and (}3T?|) we obtain the decay rate of the measured system 

R = E ~W(iO fa) = / G(co) u P(u) i>f du;, (43) 

where 

= E G ( w W ( 44 ) 

a 

The function P(u)ij is related to the measurement-induced broadening of the spectral line 
111 14,^^. For example, when instantaneous ideal measurements are performed at time 
intervals r, we can take Fi a j a '(t) = @(t — t), where Q(t) is the unit step function. Then 
from Eq. (|42|) we get 



Piool 



2 sin 2 \ \t{u — uJif] 



hf TIT (UJ - UJ if ) 2 

We have that the width of the function P(u)ij increases when the duration of the measure- 
ment t decreases. 

The equation (|3*U|) represents a universal result: the decay rate of the frequently measured 
decaying system is determined by the overlap of the reservoir coupling spectrum G(u) and 
the measurement-modified level width P(u)ij. 

Depending on the reservoir spectrum G(u) and the frequency of the measurements 1/r 
the inhibition or acceleration of the decay can be obtained. If the frequency of measurements 
is small and, consequently, the measurement-induced broadening of the spectral line is much 
smaller than the width of the reservoir coupling spectrum, the decay rate equals the decay 
rate of the unmeasured system, given by the Fermi's Golden Rule. In the intermediate region, 
when the width of the spectral line is rather small compared with the distance between u>if 
and the nearest maximum in the reservoir spectrum, the decay rate grows with increase of 
the frequency of the measurements. This results in the anti-Zeno effect. 

If the width of the spectral line is much greater compared both with the width of the 
reservoir spectrum and the distance between u^f and the centrum of the reservoir spectrum, 
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the decay rate decreases when the frequency of measurements increases. This results in the 
quantum Zeno effect. 



VII. CONCLUSIONS 

We analyze the quantum Zeno and quantum anti-Zeno effects without using any particular 
model of the measurement. The general expression (J24j) for the jump probability during the 
measurement is derived. The main assumptions, used in the derivation of Eq. (j24j) . are 
assumptions that the quantum measurement is non- demolition measurement (Eq. (JSj)) and 
the Markovian approximation for the quantum dynamics is valid (Eq. Q). We have shown 
that Eq. ()24|) is also suitable for the description of the pulsed measurements, when there are 
intervals of the measurement-free evolution between successive measurements (Eqs. ()30|) - 
(j2IIJ)). When the operator V inducing the jumps from one state to another does not depend 
on time Eq. (jSHJ), which is of the form obtained by Kofman and Kurizki is derived as 
a special case. 



[1] J. von Neumann, Mathematical Foundations of Quantum Mechanics (Princeton University 
Press, Princeton, 1955). 

[2] L. A. Khalfin, Zh. Eksp. Theor. Fiz. 33, 1371 (1958), [L. A. Khalfm, Sov. Phys.JETP 6, 1503 
(1958)]. 

[3] S. Wilkinson, C. Bharucha, M. Fischer, K. Madison, P. Morrow, Q. Niu, B. Sundaram, and 

M. Raizen, Nature (London) 387, 575 (1997). 
[4] B. Misra and E. C. G. Sudarshan, J. Math. Phys. 18, 756 (1977). 

[5] W. M. Itano, D. J. Heinzen, J. J. Bollinger, and D. J. Wineland, Phys. Rev. A 41, 2295 
(1990). 

[6] T. Petrosky, S. Tasaki, and I. Prigogine, Phys. Lett. A 151, 109 (1990). 
[7] V. Frerichs and A. Schenzle, Phys. Rev. A 44, 1962 (1991). 
[8] S. Pascazio and M. Namiki, Phys. Rev. A 50, 4582 (1994). 
[9] A. G. Kofman and G. Kurizki, Phys. Rev. A 54, R3750 (1996). 
[10] B. Kaulakys and V. Gontis, Phys. Rev. A 56, 1131 (1997). 



12 



[11] A. G. Kofman and G. Kurizki, Nature (London) 405, 546 (2000). 

[12] M. Lewenstein and K. Rzazewski, Phys. Rev. A 61, 022105 (2000). 

[13] P. Facchi, H. Nakazato, and S. Pascazio, Phys. Rev. Lett. 86, 2699 (2001). 

[14] J. Ruseckas and B. Kaulakys, Phys. Rev. A 63, 062103 (2001). 

[15] A. G. Kofman, G. Kurizki, and T. Opatrny, Phys. Rev. A 63, 042108 (2001). 

[16] M. C. Fischer, B. Gutierrez-Medina, and M. G. Raizen, Phys. Rev. Lett. 87, 040402 (2001). 

[17] V. B. Braginsky, Y. I. Vorontsov, and K. S. Thorne, Science 209, 547 (1980). 

[18] C. M. Caves, K. S. Thorne, R. W. P. Drever, V. D. Sandberg, and M. Zimmermann, Rev. 

Mod. Phys. 57, 341 (1980). 

[19] W. G. Unruh, Phys. Rev. D 19, 2888 (1979). 

[20] V. B. Braginsky and F. Y. Khalili, Rev. Mod. Phys. 68, 1 (1996). 

[21] J. Ruseckas, Phys. Rev. A 66, 012105 (2002). 

[22] C. Balzer, R. Huesmann, W. Neuhauser, and P. E. Toschek, Opt. Commun. 180, 115 (2000). 

[23] J. Ruseckas, Phys. Lett. A 291, 185 (2001). 



13 



